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ABSTRACT 



The Coppens-Sanders perturbation solution for the one- dimensional 
non-linear acoustic wave equation with dissipative term describing 
the viscous and thermal energy losses encountered in a rigid walled, 
closed tube with large length- to-diameter ratio was extended to include 
sixth order terms. The solution was then investigated to determine 
the region of validity. Computer programs were written to evaluate 
and graph the resulting waveforms. Available experimental results 
were compared with the theoretical predictions and good correlation 
was found to exist in the region of low Mach numbers. This agreement 
was found to gradually deteriorate as the Mach number was increased. 

A Fourier synthesis approach is also presented and the leading terms 



of the first ten harmonics are derived. 
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\ = 
= 



(peak) amplitude of voltage output of accelerometer 

(peak) amplitude of voltage output of microphone 

V. /A = accelerometer sensitivity 
A o ^ 



\ = 



Vj^/P = microphone sensitivity 



10 



TABLE OF CONTENTS 



Section Page 



1. 


INTRODUCTION 


13 


2. 


THEORETICAL DEVELOPMENT 


19 




Perturbation Approach 






Fourier Synthesis Approach 






Cavity Response Development 




3. 


APPLICATIONS 


4i 


4. 


RESULTS AND COMPARISONS 


43 


5. 


BIBLIOGRAPHY 


52 



3 



LIST OF TABLES 



Table 




Page 


4.1 


Tabulated Theoretical and Experimental Values of 
the Harmonic Content in the Waveforms. 


47 


4,2 


Tabulated Theoretical Values of the Harmonic Phase 
Angles , 


49 


B.l 


Symbology Pertinent to the Computer Program FINAMPI 


64 


B.2 


S3mibology Pertinent to the Computer Programs QCURVES, 
QCURC, PHAMP, and PHAMPC. 


77 



5 



LIST OF ILLUSTRATIONS 



Figure Page 

A,1 Definition of Geometrical Relationships for 56 

Both the Fourier Synthesis Approach and the 
Perturbation Approach. 

C.l Pressure Waveshapes for M = 0.004, AcO ~ O. 123 

C.2 Pressure Waveshapes for M = 0.005, 

C.3 Pressure Waveshapes for M = 0.006, 3^25 

C.4 Pressure Waveshapes for M = 0.009, ^25 

C.5 Q- curves of the Fundamental, M = 0.004 3^27 

C .6 Q- curves of the Second Harmonic, M = 0.004 128 

C.7 Q-curves of the Third Harmonic, M = 0.004 129 

C .8 Q-curves of the Fourth Hamonic, M = 0.004 130 

C.9 Q-curves of the Fifth Harmonic, M = 0.004 131 

C.IO Q-curves of the Sixth Harmonic, M = 0.004 132 

C.ll Q-curves of the Fundamental, M = 0.005 3^33 

C.12 Q-curves of the Second Harmonic, M = 0.005 

C.13 Q-curves of the Third Harmonic, M = 0.005 X 35 

C.14 Q-curves of the Fourth Harmonic, M = 0.005 136 

C.15 Q-curves of the Fifth Harmonic, M = 0.005 3 _ 3 ^ 

C.16 Q-curves of the Sixth Harmonic, M = 0.005 138 

C.17 Q-curves of the Fundamental, M = 0.009 I 39 

C.18 Q-curves of the Second Harmonic, M = 0.009 iZfO 

C.19 Q-curves of the Third Harmonic, M = 0.009 iZfi 

C.20 Q-curves of the Fourth Harmonic, M = 0.009 iZf 2 

C.21 Q-curves of the Fifth Harmonic, M = 0.009 1/^3 

C.22 Q-curves of the Sixth Harmonic, M = 0.009 

C.23 Phase Dependence of Various Harmonics on 3 ^Zf 5 



7 



ACKNOWLEDGEMENTS 



The generous aid and encouragement of Professors Alan B, Coppens 
and James V, Sanders is gratefully acknowledged. 



11 



1. INTRODUCTION 



The thermodynamic equation of state of a fluid for reversible 
acoustic processes can be written as 






( 1 . 1 ) 



or^ for an ideal gas, 

^ - ('''a) 



r 



( 1 . 2 ) 



where ^ is the ratio of specific heats. The continuity equation for 
one -dimensional wave propagation in Lagrangian coordinates is 






(1.3) 

The force equation in the absence of dissipative mechanisms has the 
form 



„ is = _ i/. . 



(1.4) 



It is possible by straightforward combination to obtain the equation 

- -z r 






fJl' 

(f). 



for a fluid obeying Eq. 1,1 or 

iif., . Y„/ „ 



(1.5) 



( 1 . 6 ) 



for an ideal gas. 
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If interest is restricted to acoustic processes for which i j 

Eqs. 1.5 and 1.6 may be approximated by 



(Mr 



(1.7) 



which (upon differentiation with respect to time) becomes the non- 
dissipative wave equation 



where 



(H) 



( 1 . 8 ) 



z 



(1.9) 



n^= .i- - ^ 

A more general form of Eq. 1.4 is 









( 1 . 10 ) 



where is an operator which generates all forces, other than that 
arising from the gradient of the pressure field, which are active in 
the system of interest [4]. If Eq. 1.10 is used instead of Eq. 1.4, 
Eq. 1.7 becomes 






( 1 . 11 ) 



In the Navier-Stokes case 













( 1 . 12 ) 



where and are the shear and bulk viscosity coefficients, respectively. 



An important contribution to the solution of Eq. 1.11 for the 
case of Eq. 1.12 was made by Fay [5]. He investigated the changes 
in the steepness of the wavefront in a periodic, finite-amplitude 
plane wave of infinite extent in a viscous medium. By means of 
Fourier analysis techniques he was able to show that the nonlinearity 
of the pressure density relationship results in accumulating distortion 
of the waveform with shifts of energy from the lower to the higher 
frequency components of the wave. Because the classical Navier- Stokes 
absorption coefficient is proportional to the square of the frequency, 
the higher frequency components are attenuated more rapidly. Fay*s 
solution is valid in a region far from the source where the rates of 
harmonic growth due to nonlinearity and the rates of attenuation due 
to viscosity are tending to balance leading to a waveform of relatively 
stable shape over distances corresponding to many wavelengths. 

Fox and Wallace [l] investigated the same problem as Fay using a 
graphical analysis technique. Their result is essentially the same 
as Fay*s. Keck and Beyer [6] developed a perturbation analysis technique 
for periodic plane progressive waves of infinite extent in a viscous 
medium with which they generate terms through sixth order. The wave 
equation is the same as that treated by Fay, but the solution is valid 
only near the source. 

Weston [20] has presented a linearized wave equation for the 
propagation of monofrequency sound in tubes. He assumes that the wave- 
fronts are planar except near the walls. The viscous and heat-conduction 
losses in the boundary layer are assumed to be the dominant dissipation 
mechanisms in the fluid. Except in this boundary layer, the wavefronts 
are independent of a radial coordinate. This suggests that the boundary 
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layer loss can be replaced by an equivalent absorptive process active 
throughout the volume of the cavity [4]. More will be said about this 
later. 

Saenger and Hudson [ 14] present a simple description of periodic 
shocks at resonance using a two part solution. Their solution is 
based on the assumption of the applicability of the linearized acoustic 
equations everywhere except in the region of the shock where the Rankine- 
Hugoniot shock conditions are assumed. They consider the case in which 
the piston oscillates at the fundamental frequency. Compressive viscosity 
is ignored but the effects of shear viscosity and heat conduction in the 
boundary layer are considered. Their solution remains finite only because 
of these considerations. 

Betchov [13] also postulates the existence of the shock wave and 
constructs a solution at resonance based on a continuous and a discontinuous 
part. For an inviscid fluid the amplitude at resonance is found to be 
finite and to be determined by nonlinear effects. The effect of wall 
friction is discussed and it is suggested that this could modify the 
solution significantly . 

Weiss [7] has applied finite-difference techniques to the nonlinear 
inviscid acoustic equations and has shown that repeated reflections at 
the rigid boundaries tend to promote the development of a discontinuity 
in the velocity profile. The absence of any dissipative mechanisms, 
however, eliminates the possibility of any steady-state standing wave 
patterns . 

While most investigations have been confined to traveling waves, 
a few have considered standing waves. 
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Solutions of the nonlinear wave equation have been obtained for 
the case of finite-amplitude standing waves in a closed tube, where 
the tube is driven at one end by a piston source. This work was done 
by Keller [2l] using the Lagrangian formulation and assigning an 
unrealistic value of -1 to the adiabatic exponent Y . Due to neglect 
of any dissipative mechanisms, his expression for the particle velocity 
becomes infinite if the piston frequency equals a natural frequency 
of the tube. 

A detailed theoretical analysis based on the nonlinear acoustic 
equations including the effects of compressive viscosity and shear 
viscosity in the boundary layer has been done by Chester [8]. He 
successfully predicted asymmetries in the finite-amplitude standing 
waves. Unfortunately his results are not easily compared to experi- 
mental results. 

Coppens and Sanders [4] developed an extension of the Keck-Beyer 
perturbation approach wherein wall losses as well as bulk losses are 
included. Their solution is applied to the case of finite-amplitudes 
standing waves in rigid walled cavities. The results yield information 
concerning the amplitudes and phases of the Fourier components of the 
waveform, and indicate the im.portance of the type of absorptive process 
on the resulting waveform. 

The purpose of this thesis is to extend the Coppens-Sanders 
perturbation approach. In part this is accomplished by a direct 
application of their method to obtain sixth-order terms. Information 
concerning the response of the cylindrical cavity to the sinusoidal 
excitation is extracted from this solution by some algebraic manipulations. 
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A Fourier synthesis approach to the problem Is presented and from it 
the leading terms of the first ten harmonics are derived. Information 
similar to that obtained from the perturbation approach is acquired 
and both results are compared to the experimental results obtained 
by Beech [19], 



IS 



2. THEORETICAL DEVELOPMENT 



In this section both a perturbation approach and a Fourier- synthesis 
approach to the solution of the one- dimensional non-linear, dissipative 
acoustic wave equation are presented. These approaches are formulated 
for a rigid-walled, closed, cylindrical cavity with large length-to- 
diameter ratio. The solutions obtained are for finite-amplitude standing 
waves. Manipulations are also presented which enable each approach to 
yield resonance response information which is convenient for comparison 
with experimental results. 



Perturbation Approach 

Assume a perturbation series such that the n-th 

term is of order n-1 in the Mach number. Substituting this into Eq, 1.8 
and collecting terms of equal order yields the set of equations 

= ^ 2 $^ ( 2 - 1 ) 



where 





( 2 . 2 ) 



If Uy) - j where the summation over j is understood to 

encompass only those frequencies included in each (My), Eq. 2.1 becomes 




7 





(2.3) 
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In order to utilize this formulation, it is first necessary 
to obtain the form of Dj. For the case of finite amplitude sound 
in a duct, Weston [20] has presented a wave equation valid for the 



propagation of monofrequency sound, Weston’s equation may be written 
in Eulerian form as 



where 



\jlu + cT U-i) 






= o 






.'So' 



(2.4) 



(2.5) 



and 

3 == (uniform) cross-sectional area 

G[ = perimeter of area S 

Cl) = (angular) frequency of the propagating waveform 
^ = kinematic shear viscosity of the fluid 

i) = thermometric conductivity of the walls 

1 . fTT 

This may be rewritten as 

- r 

When dissipation processes are included, it is usually assumed, explicitly 
or tacitly, that they are weak. Under this assumption, and for small 
Mach number, the derivatives ( ) t and (3/a X. )t may be interchanged. 

This also holds for the pair ( ^/k ) a and ( )x. Thus Eq. 2.3 must be 

solved with 



I— 1 2 /cT \ ^ ^ r « 



- o 



(2.6) 



hj 



_ j-/ 

d dCL 



1 



(2,7) 



where represents the value of cT for a frequency . 
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The final result is the approximate wave equation 

-i- 



= /, 



( 2 . 8 ) 



Eq. 2.8 was solved [4] for a uniform tube of length L excited 
at end a = 0 by a piston driven with constant acceleration at a 
frequency cJ , and terminated at a = L by a rigid cap. In this 
solution, the boundary conditions were given by 

a„(o^i) = - &,) 



where "0 will be defined later. 

The lowest order solution of Eq. 2.8 is 



a, - a,, = 






(^n) 



where 



/y - 

Un - T— ~ e 
ACd 



t{cot~e,)e - € 



ikL L 

e ~ - a 

and k = k -ic^ is the complex propagation constant. 

The solution for may be rewritten as 

^// ^ ^ (jL-^ clJJI 

Cc^ “J ,/luu j 



where 



U, =(^)(JA4U^ 



(2.9) 



( 2 . 10 ) 



( 2 . 11 ) 



( 2 . 12 ) 



(2.13) 
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and 












(2.14) 



Eqs. 2.12 and 2.13 will be rewritten in terms of (ji)^ , the 
resonant frequency which maximizes the pressure amplitude , 

associated with the rigid end of the cavity. 

lYl'it _ 

" “ L '>^7C 
cO - a) u -t A u) 






(2.15) 



The integer m represents the normal mode most strongly excited by the 
input frequency. For and 'O’f are 

2 Ao 






y^xSt cot' 



C C"3^ 



^ C7 COf^ J 



(2.16) 



(2.17) 



From this point on, higher-order solutions are obtained by first 
making the approximation 



u, 






(2.18) 



which is Eq. 2.12 with terms of omitted. The procedure (presented 

in Appendix A), is then merely an iteration process. 

In view of later discussions it is advantageous to point out that 
Eq. 2.18 may be rewritten in terms of the associated acoustic pressure 
by merely making the substitutions 

(2.19) 






Cc<3^ 
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This is equivalent to ignoring terms of order yk in Eq. 1.10, and 
converts the calculated variable from velocity to pressure, a more 
commonly measured quantity. 

In terms of the pressure, the solution presented by Coppens and 
Sanders is 



'i( ^f( 0, ( 2 . : 

+i( ^ •? 4 /j 4 I 

^ ('y '/ ^ • j 

i OA/A] + ■ ■ ■} 








■f • » • 
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where 



( 2 . 21 ) 



urar 

fitf ^ l^(Jyi~J'i T* ^ ] (2.22) 

By, = ' eot"l ( 2 - 23 ) 

M = ^“/e^ 



and M is the peak Mach number of the first order solution. 

Using the same scheme, the sixth order terms were derived. 
The solutioi) for which is to be added to Eq. 2.20 is 




, V /A7> 

Cc-v 2^(l-o~){ 



I(2,2JM>Z 



(2.25) 



-t'(2Z2,S,4)^^ + (2,2,2, SJ), 

(0 1 2 2,3), 4(1, 1, 2,2, 2),. 

3 ( /jh2,-2,2)2^ 

- i (^,2-2,3,4), y(2,2,2,-2,4)z 
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- cc^ ^ K(i-(^)i 2 f /ft - 1 



^j(2z,s,4S), 



4(z,2,3,U)_^f(^l22j4h 

^i('yh¥,‘^)* *i(2Z,¥/fh 



^[(IZ2241 




I 

-J 






*2(¥,¥,'-h 

^i(¥¥kh]4-"j 
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The entire perturbation solution herein developed, which is the 
sum of Eqs. 2.20 and 2.25, will be discussed in some detail after 
the Fourier synthesis approach has been presented. 



Fourier Synthesis 

Beginning with the approximate wave equation Eq. 2.8 and assuming 
the total solution is composed of the known classical solution 
satisfying 




^ ic 




defined by 



= o 



(2.26) 



/c ^ 



(2.27) 



then 




^ _ /• 

-1- tt' ^ 



(2.28) 



Define the phase velocity 




to be - Co^ ( h dj ) 



and Eq. 2.28 may 



be written 



[cJ C St ^60^) 

”” Z \ ^ 



(2.29) 



Assume ^ may be written as a Fourier series 

Xy. 0<3^(y\0St+ (/) 



(2.30) 
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and 



^ = Xc Co^ Co^ cot 

and let 

(^f=£A(<^) co^y<cot fO, ) 

1\ -c 

where X>) and />) will be obtained later. Combination of Eqs. 
through 2.32 yields 

L(cV^~ ^ 0^ ”) 

- cfvj X>T 'y\K(L-&S) Ji^ (^ost + 



We impose the conditions that 

X(V = ^->1 Co^ 'hk^Cl-a.) 

and 

/jy^'fo.) o<: ^y\ K y\!^(L-OL) 

so that Eq, 2,33 becomes 

= ^ 3>i ) 

Introduction of the phase angles of Eq. 2.17 allows Eq. 2 
written 

-Y^ f s ^ Co^ 2L -f-/^) 



(2.31) 

(2.32) 
2.29 

(2.33) 

(2.34) 

(2.35) 
•35 to be 

(2.36) 
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In order to simplify the manipulations to obtain and 

define 

X,Co^k.(i-CL ) C<yi(oA-t4‘i) - vrcoi K(L-0^) Cen(c^i'-f- °) 

C<y-v k(L~o^) Co^CO't 



and 



Thus 






Xh K(L~Cl) Co^ (^U)t ) 



yt-l 



and 



Cc^ 

<x> 

==■ y~^ \-j Ccru (huft -f ~ 9^^' ) 



r 

Combining Eqs . 2.40 and 2.36 gives 



•4^ 6^ 

Cc4 



Xv» - B-in) 



CO 

X■^y^h■^■j C<n(ho)i: ^ 



■= / 



-h-l 



^ -i- 1 






(2.37) 



(2.38) 



(2.39) 



(2.40) 



(2.41) 
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It is convenient to shift to complex notation at this time. Let 



X, - X,° 4>° 


(2.42) 


SO that 


(2.43) 


and 

Ay>e = X>, e , 1 


(2.44) 


Then Eq. 2.41 becomes 

y 0 c ■&T 1 -'^z) 

(-ri C<yi^-&y\ 

r' /-/ 

If the equivalent pressure equation is desired, make the 
substitution 


(2.45) 


X — n = - 


(2.46) 



where M is the amplitude of the Mach number for the particular harmonic 
Define 



- r/2tly, 


(2.47) 


SO that Eq. 2.52 becomes 

JV -Vz) 

■St, 

= 23";^” '/ ^ ^ 
2T 

r' 


(2.48) 
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Examination of Eq. 2.48 and Eqs . 2.20 and 2.25 shows that the two 
methods are equivalent since expansion of the term 

i ^ Mj 

yields the leading terms of each harmonic derived in the perturbation 
scheme, while the other portion of Eq. 2.48 contains the correction 
terms . 

Recall that in both the perturbation approach and the Fourier 

synthesis approach the leading term of the fundamental was calculated 

to be , This is of order one in the Mach number. 

It is also a known experimental fact that the fundamental remains 

essentially unchanged with respect to Mach number. This suggests that 

the phase angles associated with the corrective terms might cause the 

corrective terms to sum to approximately zero. It then becomes 

conceivable that a solution composed of only leading terms could 

yield significant results. In order to investigate this possibility, 

/ 

Eq. 2.48 was used to obtain the leading terms of the first ten harmonics. 
The result is 

+ J (!^) (Z]A 

-/• j Cc-u 4 {(0[wu 
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* i ( imwJf 

-^/i(^.w,¥^‘^)f]j 



// (’«. j(^/f(2j¥:ri7/J/ol 



fo 



fc 



(0 



\ 

*2 A ^ Mi¥.i^zn'o), 

+i (2234X?fiVo), 

AA//,^A/A a umv.v^i 



/6 



tC 



/C 



fC 
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+^ (22234589/0), 



i 



i ( l^,VA^7,V°),c ^ j (22 3 3 44 89, /o), 
± 



/O 



* K '^^7777A^'°),> 4^(22uHe9/o),^ 

^77/°)'- * j(22245i79/c),„ 
"^4 ^(^22 3^^ ^7^9/cJ^ 

^/i ^777 ^ JHo " ^(2223U^f/o)^ 

^ s^^^^7V3'‘'),o ^ i(777,¥,V/i. 
itU'(W7AV‘‘),c 2 
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^ is 7?.W7M'4,, * i(223s,n^7/4/o 
4 i (9.77773/4,0-^ )(99A7Ho 

■^8 (77/7, 9/4,0 ^ 2^(99 779/4o 
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7<r 



7<? 
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Vci 



7<l 



fc 



/c 



'*'8 

* ^ <' ^ /if' 7l^^-fA78/‘), 

* i ( ^A^.VAA,‘SmI^ a i (^AAAAps/c), 
A \ ii . WAFA ^ A^k + i ^^ AA . WAA /^). 
/ s(^AAAAAAfj% f i (2. >.},■>, i,V‘J. 

^jz ^/i 

>'/ < A , ^ A 4 (^ AAAAkA^Ai 

i { 2 AAJ . WA /“}, 

\ 

+ i QaM j^q{WAA¥A'^X 
-^i(^,VA//A, 1'“),, * j (272}As;i 7/oJ, 

^4 (^AA/A/AVXj, -^^(222S44^7/cI 

+ji (27,72 3Uj/ol^ s (277 IpA ^ 'X, 
(2,7,7 mV^k ^ fi(277734/,7/c), 



/« 
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^ji (^.WM^A< + i(mi^,w^),> 

* 4- ( ■^/i(WAAA 
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/c 



'/c 



7c 



where (a,b)g has been given in Eq. 2.21. 
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Cavity Resonance Response 



Investigation of the behavior near resonance of each approach 



leads to the desirability of obtaining the response of the cavity 



for each harmonic as the fundamental is s\\;ept through resonance. 

It is obvious from the form of Eqs. 2.27, 2.32, and 2.55 that both 
approaches yield results which may be written as sums of terms 
having like factors of K(L-a) andcO t . It is also evident that each 
of these harmonics is heavily dependent upon M, the Mach number, and 



The algebraic process necessary to bring each of the two solutions 
to a form which conveniently yields the desired information will be 
derived in a general form and then applied to the perturbation 
solution. It is then easily extended to the Fourier synthesis solution. 
We begin with the relationship 

^ Qc 4c) — ^ -^ (^ ) , (2.50) 

c' = / 

At t = 0 Eq. 2.50 becomes 
L 

The derivative of Eq. 2.50 evaluated at t = 0 yields 

'yi 

~ (2.52) 

c-( 

which when combined with Eq. 2.51, yields 







(2.53) 



If Eq. 2.51 is divided by Eq. 2.52 the result is 




(2.54) 
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The sum of Eqs . 2.20 and 2.25 may be written as 




(2.55) 



for each harmonic, where j indicates the harmonic number and n 
indicates the number of terms in that harmonic. It is evident 
that the amplitude of each harmonic may be obtained by applying 
Eq. 2.53. 

In applying Eq. 2.53 to the perturbation approach the magnitude 
of each harmonic was normalized to approximately one by the maximum 
magnitude of a factor chosen from the leading term of that harmonic. 
This normalization was done to yield a set of curves whose maxima are 
of order unity. For the fundamental this normalization factor is one. 
The method of choosing the normalization factor &r higher harmonics 
will be demonstrated for the second harmonic. The leading term of 
the second harmonic is 




(2.56) 



Now 




(2.57) 



so that 




(2.58) 



and the leading term of the second harmonic may be written 




(2.59) 
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Since 



Eq. 2.59 may be v?ritten 

- 4 ^ 

("o 41 



H has as its maximum 
n 




(2.60) 



(2.61) 



(2.62) 



and the normalization factor for the second harmonic is then chosen as 



- 4 ^- 

27*)7t^f,C(, 4/ OT- 

It should be noted that these normalization factors are not the maximum 
amplitudes of the various harmonics, neither are they the maximum 
amplitudes of the leading terms of each harmonic; rather, they are 
factors chosen to yield maxima of order unity. 

These factors are given the symbol MAXA^ where n refers to the n*-th 
harmonic so that 



= 





( 2 . 64 ) 



By combining Eqs . 2.53 and 2.64 it is possible to define as 
normalized response of the n-th harmonic. It should be noted that A 
defined in Eq. 2.53 is a function of the Mach number and . 

Thus 







(2.65) 



the 



where n is the harmonic number. 
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3. APPLICATIONS 



Associated with each of the solutions formulated in the previous 
section are three computer programs. These programs are concerned 
with the pressure waveform, the Q-curves, and the phase information. 

The perturbation solution has the programs FINAMPI, QCURVES, and PHAMP 
connected with it while FINAMPIZ, QCURC and PHAMPC are related to the 
Fourier synthesis . 

The first attempted utilization of the perturbation solution 
consisted of writing a computer program involving all terms through 
sixth order. This program, entitled FINAMPI, was designed to compute 
and graph both velocity and pressure waveforms. FINAMPI was first used 
to graph the waveforms contained in Ref. [4] and has been used since 
in determining for what Mach numbers the theory begins to deviate 
significantly from the experimental results. The graphical output 
from the computer is contained in Figs. C.l through C.4. 

The program designed to compute and graph the cavity resonance 
response for the perturbation solution is entitled QCURVES. This program 
is a direct application of Eq. 2.65 and has been used to obtain the 
curves contained in Figs. C.5 through C.22. 

PHAMP, the third program associated with the perturbation solution, 
was designed to compute the phase angles of each harmonic as the 
fundamental was swept through resonance. It is a direct application 
of Eq. 2.54 and resulted in the data contained in Table 4.2. 

In an attempt both to determine the importance of the corrective 
factors to each harmonic and to test the validity of the Fourier synthesis, 
computer programs involving only the leading terms of the various harmonics 
were written. Unfortunately a program involving all ten harmonics could 
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not be brought to fruition because of time limitations. A second 
program, FINAMPIZ, which uses the leading terms of the first six 
harmonics has been used to graph the pressure waveforms also contained 
in Figs. C.l through C.4. 

Eqs. 2.65 and 2.54 were also applied to the Fourier synthesis 
yielding the programs QCURC and PHAMPC. QCURC computes and graphs 
the Q-curves for each of the harmonics while PHAMPC computes and graphs 
the phase dependence of the individual harmonics on A cO . 

It should be noted that of the two programs concerned with phase 
dependence, only PHAMPC graphs its output. This was done because the 
Q-curve programs indicated the Fourier synthesis approach produced 
curves which were in much better agreement with the experimental results 
than did the perturbation approach. The graphical output of PHAMPC is 
contained in Fig. C.23. 

In this section only a superficial description of each program 
was given. A detailed description of all programs is contained in 
Appendix B, and a more complete discussion of results follows in 
Section 4. 
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4, RESULTS AND COMPARISONS 



In this section the theoretical predictions obtained from the 
applications described in Section 3 will be compared to the experimental 
results obtained by Beech [l9]. 

Before proceeding to a discussion of these results it is necessary 
to attempt some qualitative evaluation of the apparent properties of 
the perturbation solution. The amplitudes of the (n + l)th order terms 
tend to behave as [4]. 



01 



iflt i = 

i’Z / 7'^ 



(4.1) 



It is clear that this expression can diverge for increasing n. This 
divergence, which occurs because of the dependence , would 

not occur if since the left hand side of Eq. 4,1 would then 

be bounded by ^ series would converge for /^ < 

The Mach number which satisfies the expression and would 

therefore appear to be related to shock formation is M = 0.02. This is 
within the same order of magnitude as the experimentally measured Mach 
number at shock formation, M = 0.01. 

The divergence may be remedied if it is recalled that the bulk 
absorptive processes have been neglected up to this point. If the 
viscous terms are retained and assumed to be additive to the wall losses 
then 




(4.2) 



where 




(4.3) 
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The analysis and results are the same as before with the H*s and 0's 



redefined as 



/ -y\cd^ 









(4.4) 

(4.5) 



With these modifications the amplitudes of the (n + l)th order terms 
now behave, for sufficiently large harmonic, as 

pt f 

which tends to zero with increasing n. 

The modifications introduced by Eqs . 4.2 and 4.3 are important 

in the perturbation expansion for those harmonics whose indices 

J> / . 

satisfy the inequality ^ /cO^ . For the problem dealt with 

in this research the theory developed without the bulk absorptive 
terms is adequate as long as attention is restricted to frequencies 
below about the thousandth overtone. 

In comparing the theoretical and the experimental results it is 
convenient to divide the discussion into three areas: the low Mach 

number region, the intermediate Mach number region and the high Mach 
number region. 

In the low Mach number region, below M^ 0,006, it is obvious 
from. Figs. C.l and C.2 that theory and experiment are in very close 
agreement. In these figures of pressure waveforms the large drawing 
contains the waveforms predicted from the computer programs and the 
smaller inset is a drawing from an oscilloscope photograph of the 
waveform observed experimentally. The dashed line in both drawings 
is a sine wave which is included for reference purposes. One of the 
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main reasons for the good agreement in this region is due to the 



It should be noted here that in all the figures of pressure 
waveforms, the major difference between the FINAMPI and FINAMPIZ 
predictions is in the phase dependence. Recall that the development 
of the reasoning behind the FINAMPIZ program rested on the possibility 
that corrective terms to the generated harmonics could be ignored. 
These curves lend support to that suggestion in the low Mach number 
region. 

The low Mach number region being discussed also includes the 
first two sets of Q-curves, Figs. C.5 through C.16. It is evident 
in both sets that the agreement of the experimental results with the 
predictions of the Fourier synthesis approach is excellent in the 
lower harmonics and deviates only slightly in the higher harmonics. 
When the agreement does deviate, the experimental results in general 
fall between the two theoretical predictions. This tends to indicate 
that the two solutions have the exact solution bracketed. It is 
evident in the QCURVES predictions for M = 0.005 that the third 
l^armonic curve is no longer a smooth bell-shaped curve. The reason 
for this is the growth in the third harmonic of the correction terms 
arising from the fifth-order solution. For instance, the ratio of 
the approximate magnitude of the first corrective term to the leading 
term in the third harmonic is about 0.5 at M = 0.005. This is no 
longer insignificant so that the predictions based on the perturbation 
approach begin to show unrealistic fluctuations. 



fact that the 




terms are as yet relatively small. 
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Theoretically predicted amplitude information is obtained from 
the two programs, PHAMP and PHAMPC, The amplitude results from these 
programs together with the experimental results are contained in 
Table 4.1. It should be noted that Beech's uncertainty of + 0.1 Hz 
in frequency results in an uncertainty of + 0.628 radians in Figs. C.5 
through C.23. A comparison of the uncertainty in amplitude (based 
on the uncertainty in AtD) with the observed discrepancy between 
measurement and theory for M = 0.004 anddo0= 0 yields 



n uncertainty in 

harmonic content 



discrepancy 
between PHAMPC 
and experiment 



2 



0.5 



0.4 



3 



0.3 



0.17 



4 



0.15 



0.12 



5 



0.07 



0.08 



6 



0.03 



O'. 07 



It is evident that the uncertainties in harmonic content are consistent 
with the discrepancies in harmonic content noted between PHAMPC and 
experiment . 

As was stated in Section 3, PHAMPC also predicts the phase 
dependence of the various harmonics . The results of this program 
are contained in Table 4.2 and Fig. C.23. One thing should be noted. 
The sets of phase angles predicted by PHAMPC are identical for all 
Mach numbers, which is not surprising if the reader recalls that the 
normalization process destroyed the Mach-number dependence of the 
leading terms of each harmonic. 
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Tabulated Theoretical and Experimental Values of the Harmonic Content in the Waveforms = 



One other program exists for the investigation of phase and 
amplitude characteristics. The computer program FOUANAL is being 
introduced here since it was not derived from either theoretical 
approach. FOUANAL takes data from an oscilloscope photograph and 
computes the Fourier coefficients of the various terms of a Fourier 
coefficients of the various terms of a Fourier expansion of the 
waveform. The results from this program are also contained in 
Tables 4.1 and 4.2. 

It is clear that the results from FOUANAL are fairly accurate 
in the lower harmonics (first, second and third) but tend to become 
much worse in the higher harmonics (fourth, fifth and sixth). Judging 
from the amplitude results it does not seem likely that FOUANAL would 
yield accurate results in the phase angles of these higher harmonics. 
Examination of Table 4.2 verifies this supposition. These erratic 
results in the FOUANAL calculations may possibly be explained in the 
following manner. The input data for the program come from extracting 
sixty- four points from a standard oscilloscope photograph. This was 
done by hand and was subject to many errors. It is possible that the 
data were not given to sufficient accuracy to enable the computer to 
predict accurately past the third harmonic. The most feasible suggestion 
for improvement in this respect is to use a digital counting procedure 
to obtain the data. It should be noted that the data in Fig. C.23 
are the FOUANAL results for M = 0.004. Clearly these points are 
extremely random and the uncertainty in the experimental cannot 

possibly account for the entire discrepancy. 

Shifting now to the medium Mach number region (0 . 006 0 . 009) 
we see from Fig. C.3, for M = 0.007, that although the general agreement 
is still good, deviations from the previous good agreement are becoming 
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Tabulated Theoretical Values of the Harmonic Phase Angles 



evident in the computer predictions. Unfortunately experimental 
data for Q-curve and phase angle comparisons were not obtained 
but it is not unreasonable to assume that these comparisons would 
show essentially the same result as stated for the pressure waveform. 

Moving into the region of high Mach numbers, M = 0.009, we find 
that the pressure waveforms predicted by each program are generally 
invalid. This is nearing the experimentally observed Mach number for 
shock formation (M = 0,01), and since the theory is predicted on small 
Mach number, pre-shock conditions it is not surprising to see these 
deviations , 

Essentially the same results are observed in the Q-curve comparisons. 
Figs. C.17 through C.22. The reason for the experimental amplitudes 
falling off so rapidly is at present unknown. Tables A.l and 4.2 
again reflect the serious discrepancies in the high Mach number region. 

It is now possible to make some general comments concerning these 
two approaches. Both the perturbation approach and the Fourier synthesis 
approach have adequately described the effects of finite-amplitude standing 
waves in rigid-walled cavities provided the Mach number is relatively low. 
This precludes that area in the pre-shock region where significant 
distortion exists. It has been demonstrated that in the region where 
the approaches are fairly accurate, Q-curve information can be easily 
obtained. This information is convenient for comparisons with experi- 
mental results. It is also possible to obtain some information on the 
phase angles of the harmonics. 
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It seems plausible now to suggest the following: Since the 

Fourier synthesis approach gave the most accurate information, the 
computer program involving all ten harmonics should be reviewed in 
the hope that a usable program will result. This program could be 
used to investigate the effects of the higher order harmonics on 
the existing solution and possibly give an indication of whether 
extending the present developments past six harmonics is beneficial. 
Finally, after some exposure to basic communications theory, it 
appears that a good grounding in this subject is mandatory if the 
present theoretical line is to be pursued further. It would not be 
surprising to find that a more exact solution may be obtained as a 
result of the techniques used in communications theory. 
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APPENDIX A 



Coppens-Sanders Iteration Scheme 

From Section 2, the approximate wave equation for f inite»amplitude 
standing waves in a rigid walled tube is 






r “'7 



2 yi-l 



dh-i . V 
" "SET 9-hj 

7 

For convenience the right hand side of Eq. A.l will be written 



(A.l) 






Assume 



7 



/ 



(A. 2) 



-V (A. 3) 

where the and are complex, and the three dots in the 

exponential term represent phase terms to be defined. Work is done 
in one n at a time so that the n subscript may be dropped. 

Expansion of the left hand side of Eq . A.l yields 






(A. 4) 



where 



(/-f; 



was previously defined to be 






Thus Eq. A.l becomes 

. z -“2 






(A. 5) 



for a given order n. If Eq, A, 5 is generalized in time it becomes 

iA- J. ^ / 

(A. 6) 
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Substitution of Eq. A. 3 into Eq , A. 6 yields 



Now each ^ will be of the form 

^ ^ Cou ciJ'^ ^ ^ ) 



Substitution yields 






(A. 7) 



(A. 8) 



(A. 9) 



Additional substitution of Eq. A, 9 and Eq. A. 3 into Eq. A. 7 gives 

- /: / ^ 



uA 



1 







aJ 



(A. 10) 

Proper manipulation of the bracketed term in Eq. A. 10 will yield the 
phase relationships mentioned previously but left undefined. 

First recall that 



CO =" -f Ao) 

(A. 11) 

-2>rr __ _ y_ ^ 

c ^ M 

where 

XT = ^ 

= ^7' 

and notice that 

^ ^ A'W'^ ^ (A-12) 

and 

^ = tT (A.13) 
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(A. 14) 



Substitution into Eq. A. 10 yields 






cOf. 




It is now possible to write down the now U as 












Define the following quantities; 




and 









Pictorially these definitions may be represented as in Fig. 




Pictorial representation of the quantities 
defined in Eqs. A. 17 through A, 19. 

Figure A . 1 



(A. 15) 

(A. 16) 

(A. 17) 
(A. 18) 

(A. 19) 

A.l. 
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With these definitions, Eq, A, 16 becomes 




and 






^7 '^7 



(TV* 



yielding 





~ e ^ ^ ^ 



(A. 20) 



(A. 21) 



(A. 22) 



The preceding was a general outline of the procedure to be 
followed in the Coppens-Sanders iteration scheme. The remainder of 
this section contains a derivation of the terms through third order. 
Assume the input 

^ (*-23) 

which was shown in Section 2 to be the calculated first-order solution. 
Then we have 









£j 

and since 
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(A. 24) 
(A, 25) 






(A. 26) 



The next step is to obtain * 

~ ^ - <5^ /*» 

c)cu 

6 ' / 



i^z ' Zr 






(A. 27) 
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Substitution of Eq. A. 26 into Eq. A. 27 yields 



- /V /(' ( /.-a. ) yduK, lij 

where M = ^"/Co 

Expansion of the trigonometric term gives 

(?<«- ^ ^ uft 



(A. 28) 



(A. 29) 



== ‘ 1 / - / 






from which the constants + 1 and - 1 may be discarded since application 
of will remove them. Then 



which implies 

A. ^ _ 



/f2t = - / /V ^ 



(A, 30) 



(A. 31) 



and from Eq. A. 20 

< 5 ^; - 

Recall that hence, 



(A. 32) 



{^j /i -^^ 2 ) (A. 33) 

which is the second order solution. 

The general equation for is 

^ TV! (A. 34) 

Substitution of Eq. A. 34 into Eq. A. 21 results in trigonometric products 
such as 

CcO- ^ K>i) Ccr< ^ J (A 35) 
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(A. 36) 



For simplicity let 

X ? ^U->. ) 

/ -f 

Then Eq. A. 35 may be written 

(Pen CrL JX 

Expansion of Eq* A. 37 gives 

C~ ^ 

- c<H. (/n-j)x y 

where terms of the form sin k!(lr a) cos j^coi: have been ignored 

in Eq. A* 38 since substitution of these terms into the perturbed wave 
equation result in amplitudes of order Vi which have been neglected 
throughout. Eq. A. 38 must always be written as (m - 1) where m ^ 1, 
since a positive frequency has always been assumed. 

Before preceding to the derivation of the third harmonic terms, 
a brief discussion might be of some use in acquainting the reader with 
higher order terms. 

In the third harmonic derivation, the user of this iteration 
scheme is first exposed to the development of more than one term in 
the solution for a particular order. In other words, the iteration 
yields both the leading term to the third harmonic plus the first 
correction term to the present solution. At the same time, the use 
of Eq. A. 38 is demonstrated for the first time. It should also be 
pointed out that the methods used in the third harmonic iteration are 
identical to those used in the higher order iterations, regardless of 
the number of terms encountered. 



(A. 37) 



(A. 38) 
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(A. 39) 



Following the step by step procedure then, 



Vj - 7 ^ ^ -h ^ 

,5»tx cj^cc c)c^ 

(T'"/ 



Remember , 

/; ^ ^ JLi^ ojt 



y^> :sr-2^ Qr^ 






and 




j. 

//z 



<f 






SO that 



Tpi -2(- ^ ^ e^ 2U- 



a 




>f* 






or using Eq. A. 38, 



— a*o±. A^(/~<^J Co^ &z)]_ 



(A. 24) 
(A. 25) 

(A. 40) 

(A. 41) 

(A. 32) 

(A. 42) 
(A. 43) 
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Using Eq. A. 8 implies 




(A. 44) 



(A. 45) 



and from Eq. A. 10, 








Z ^ 

/<4 //j 



^ c(e,i^s) 



and 

list 



_ _l -rj/M)^^ - 



(A. 46) 



(A. 47) 



Now > and when the substitutions are made, the 

final form of the third harmonic iteration becomes 




6l 



A combination of the first three orders to obtain the total 



solution thus far results in 



Notice that the second term of the third order solution is of order 
one in j, (with reference to Eq, A. 3), and thus it is termed the 
first correction term to the fundamental or first order solution. 

Rewriting Eq. A. 49 in order to achieve a form consistent with 
this last observation and the form of Eq. 2.20 results in 







-iMfJ- 



Coi^o/ ^ ^ ^fi~ j 



-f-' 




-t 






1 

//z 







/ 




(A, 50) 



/ • J 



^ 9 < 
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APPENDIX B 



Computer Programs 

An IBM- 360 digital computer v?as utilized in performing the 
calculations, and obtaining the v;ave forms predicted by the theory. 
Computer programs v?ere written for each of the basic formulations 
as well as for each of the variations invest igated . Each of these 
programs, together with a short description, is presented in this 
section. Included with the description is a list of the symbology 
peculiar to that particular program. 

B-1: Program FINAMPI 

This program has written to perform the calculations and plot 
the graphs for the perturbation solution. It was originally written 
in FORTRAN 60 language and subsequently converted when the IBM- 360 
became available. 

The program was designed so that the Mach number, resonance 



are the input parameters. The first half of the program is devoted 
to the velocity profiles and the second half to the pressure waveforms. 
The functions defined at the end of the main program are part of the 
conversion from FORTRAN 60 to FORTRAN IV, A list of the symbology 
pertinent to this program, together with the corresponding parameter 
in the original formulation follows. 



frequency, length of the tube, position 




and 
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Table B.l 



D(N) 


--<5^ 






S 


= Mach number 






WD 


= Aoj 






WR 


= (Oft 






B 


II 






F 


= L 






A 


= a 






U 


= Velocity 






P 


= Pressure 






PI 


= 'a 






WT 


II 






Q(K) 


= Sk-y2iik 






THET(N) 


= 






Vxyz 


= the expansion of 


<x,y)^ 


for the velocity equation 


Yxyz 


= the expansion of 


(x.y). 


for the pressure equation 
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B-2: Program F INAMP IZ 



This program was designed to duplicate the output of FINAMPI 
using the Fourier synthesis formulation through sixth order terms. 
Since the final results of the perturbation approach are in the same 
form as those of the Fourier synthesis approach the program FINAMPI 
was converted to FINAMPIZ by setting the corrective terms to zero. 
Because FINAMPIZ is essentially FINAMPI, Table B.l is applicable here. 
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B-3 : Program QCURVES 



This program was designed to compute and graph the normalized 
Q-curve for each of the six harmonics using the perturbation solution. 

In this program, the input parameters consist of the resonant 
frequency, the tube length, ^ and the piston acceleration. The 

symbology pertinent to this program is given in the following table. 



Tab le B . 2 
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SMALA 
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piston acceleration 
length of tube 

hypothenuse of triangle described in Fig. A.l 

unnormalized magnitude of the n-th harmonic 
at point m 

the normalizing factor defined for each 
harmonic 

defined subroutine to main program 
defined subprogram to main program 
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B-4: Program QCURC 



This program is the Q-curve program based upon the Fourier synthesis 
approach using only the first six harmonics. The curves are normalized 
in exactly the same manner as in QCURVES and no new quantities have been 
defined. 
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B-5 : Program PHAMP 

This program was designed to compute the amplitudes and associated 
phase angles for each harmonic. Essentially this program is the QCURVES 
program with only minor modifications. One new quantity has been defined 
as PHASE(I,J). This is the phase angle associated with the i-th harmonic 
and the j-th in all other respects Table B.2 is applicable. 
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THIS POKTION or THF PROGRAM COMPUTTS THE AMPLITUDE AND PHASE HE 
THE FIRST HARMONIC 
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WRITE ( ft, 210 ) 

210 FORMAT! IHO ,5X , ]HM, 1 C X , ? HE I , 1 P V , ^ h r ^ , I 7X , OHF ACT ( 1 , V ) , / ) 
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THIS PflKTION r.F THE PROGRAM CUM 0 (JTFS THE aMPLIT'.nP AND PHASE OF 
THE SIXTH HARMONIC. 
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B-6: Program PHAMPC 



This program is identical to PHAMP except that it was written for 
the Fourier synthesis approach and graphs the phase angles associated 
with each harmonic. In order to use the graphing subroutine the 
subscripts on the phase angle had to be reversed. PHASE(I,J) is now 
defined as the phase angle associated with the i-th AcO and the j-th 
harmonic. In all other respects Table B.2 is applicable. 
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B-6 : Program PHAMPC 



This program is identical to PHAMP except that it was written for 
the Fourier synthesis approach and graphs the phase angles associated 
with each harmonic. In order to use the graphing subroutine the 
subscripts on the phase angle had to be reversed. PHASE(I,J) is now 
defined as the phase angle associated with the i-th AcD and the j-th 
harmonic. In all other respects Table B.2 is applicable. 



Ill 



I 



( 



>• 








, 










• CD 


«0CD 
















UJ 


















OJO O 


UCD 




X 












OOZ UJ 


CLO 




O' 












<<Z»Ki 






•» 












XZ<*-^ 




•» 


o 












CLO 




» 


< 












U0O< 


00 OD 




X 












OUJZX 






CM 












ZQ^<Q^ 






m 












< o 






X 












xoz 


(SO •►nO 


m 


00 












UJOUJ • 




oc 


•» 








•ii* 




03MZ0 


CD«— 




X 








L- 






aoooo 


X 


X 












h*o:.juJZ 




o* 


*i"4 








UJ 




•-X<COO 


CJiOUJ 


CO 










o 




-JH-X X 


O3C00 


m m 


X 








X 




a Q^OOQ^ 


00 •>< 


u. 










r- 




XKO<< 




UL II 










•> . 




<Q-ZXX 


— -4Q. 


• o 


< 








X 




UJ 


^00 




►- 








o 




UJ3C<Xh- 


00^-^ 


m 


UJ 








•"M 




Xt>n QC< 




3 


CO 












I- XOX 


<xoo 


QC 


X 








LL 




uOh-CLh- 


CL< ► 












d. 




UJi— O 


*«X^ 


mm m 










X 




CDOLL 




o ^ 


X 


u. 






< 




3V UJO 




• • ll 


h- 


m 






X 






00 *-40 


• » 




o 










XZ^>^X 


•‘Q0< 


a o 


QC 


< 






m 




OUJ 


-^wLL 


- X < 


X 


•» 






X 




U3^-<UJ 


h-u. 




X 


X 






o 


•-If*. 


03Xh- 




- < X 


(\J 








•"M 




oujoa: 




K •“ 


» 


< 


#«* 








►— 0^ 


•>0^00 


-J < 


X 


K 


O 




o 




U-OZOO 


— 


OUJ H 




UJ 


o 


CL 


X 


.JU 


O UJ 0^ 


r^OOLL 


Q.« • UJ 




CO 


* 


X 


•"M 


Ui 


UJUJKUJ'-< 




X »co 


««» 


•k 


q: 


Ui 




oo 


ZXZXU- 




^m m ^ 


•"Ml 


QC 


X 


H 


X 






UJ-4-^ 


XOUQ^ 




X 


♦ 


•■* 


o 


u • 


^ Ct UJ 


XCO*-^ 


Q. —X 


>- 


•» 


M 


H* 


•"*4 


OLcn 


C0(/)Q. »X 


►—*00 


• OOZ •“ 




— *k 


a. 


QC 




X-4 


UJ< X>- 


*.UJ-*LL 




UJ 


•"M 


# 


a 


tai4 


<u 


O 1>0Q^ 


-*1^0. UJQ. 


— ^ X 


o 




X 


00 


X 


•> » 


0»-OLL 


'f) — OX 


(NJUJOlw 


X 


L- 


w> 






OlTN 


O0I-- U.O 


---QC 






O 








» • 


^ZUJ 




^OX-J-^ 


p> 


UJ*>» 


0>^ 




X 




OOOUJ 


^oo-«X 


XUJ* • ujin 


X 


om 


< 




o 


*m4 


XX3UJO 


UJ »oo 


J • 


(NJ 


• 




L- 


*m4 


-‘U 


<Cth-KJ3 


o »- 


V.K* • O 




-.o 


UJO 


O 


un ^>45^^ 




*-^OUJ 


jN-»z 




CD » 


Ui 


CO ^00 


OX-J-J^ 


zcoaxx<ujK< oll-41 


rsjtu^ • 


♦ -4 


O 


•«"4fi^ 


O CL<-J 


o— •^UXCDi— -^sO 






O II 


II UJ 


•.•-4 (1 


•1-1 UJ 


o:xxxcL 




«< UJ ^ 


nO-— 


<1-4 




IO--3 


Q-0<Q^X 


ooQ.— 4ZXX— iooco ir^L- 


h- ^ ^ < fO LL -* 


>-4 NK Z 




— H-Z 


< o< 


Z«-4 00O 


•»< — <UJ< 


•»UJ< • 


M II 


uj<h-uj<#- 


OOUJXZ 


UJO^^X 


-J-JU-JX OX>-XU-KXrO II 


^y^*-4 








« O • UJ 


xoox<<<<a <QC*--ocx*-x II 


aa 


X-JZ 


«QC 


•-QCZ 


XU-<ZX 


^ooOOUJUJUJUJ» • UJ0CiC0-<CiC0'-»0 II XOUJUJOQCOOQ^OO 


►•OUJ3^- 


O i^QCOCaCOC OlLL 


XLL 


XU.Q"00<QH'OOXU-DXLLO 






— 4fSJ 


9 


•"4 














O 


U) 








ITk 










•m4 


CVJ 




*m4 




00 ir^ 



oouooooou 



ouo 



112 



THIS PORTION COMPUTES AND GRAPHS THE NORMALIZED Q CURVE 



o 













•» 










X 




CO 










CVj 














O' 




•k 










* 




X 










•k 


















w 










««« 




•k 














X 










X 




o 


UJ 








X 




X 












kkifc 


•k 


> 






















< 




o 


QC 






«*ik 


1- 














X 


•k 


•k 


Z> 






fVJ 






X 










< 


X 


o 


o 






« 


Uj 




in 










X 












* 


o 














X 


UJ 


o 


Of 










<»• 














00 


•k 








X 


+ 


* 












•k 


< 


o 


a 






w* 






X 










X 


X 




UJ 






o 


fSJ 














r- 


CL 


o 


INI 






!5 


•w* 




h- 










•k 






Pm4 






W 




>- 


UJ 










k-k 


k-* 


tu 


«J 






1 






X 










X 


X 


mJ 


< 










UJ 


H 










W 


•Wk 


K 


X 






0^ 


or 


c 


X 










<r 


< 


M4 


oc 






3 


3: 


w 












CL 


a 


1- 


o 






•«• 




V. 












— t 






z 












X 










QC 


kM, 


•k 








X 


X 




>c 


X 








U 


X 




UJ 








w 


• 




o 








to 




JJ 


X 






o 


o 






< 








X 


< 


X 


1— 






3T 




w 


CO 


X 








cr 


X 


< 








# 


« 


K 


X 










•k 


< 




00 






• 


• 


a: 












X 


X 


«— «k 


X 






<\J 


Cvj-*> 


o 




X 








in 


•k 


o> o 


a 








w<^ 


i/> 


X 


w 








•kk— k 


*— 


W *k 


< 






1 


1 ^ 


•«« 


CO 


X 










X 


UJ O 


QC 








»■* w 














X 




-J 


O 






rsj 


XX 




<► 












< 


K — 














< 


X 








Of-4 


a 




O 






* 


1-0 


* 


X 










3C 


»«« 


K 


z • 






a: 


uo 


«i>% 




►- 








XX 


QC 


— * 


<o 


o 


in 


3L 


UJ«- 


< 




Lli 




X 




inw 


O 


•k w 


1—4 




nj 




oa: 


K- 


X 


X 


w 


w 




•‘UJ 


00 


^ Oi 


002 


2 




* 


1 :x 


UJ 


o 


K- 


1— 


< 


GO 


XoO 


•k 


00 00 


UJO 


O 


• 






OD 




•» 


UJ 


X 


•» 


oo< 


— 


Mk < 


l-X 


T 


o 






■«• 




GO 


X 


< 


< 


•^X 


X 


UJ X 


OQC 


Ct 




X— 


00*^0. 


X 


X 




K 


X 


w 


XCL 


W — «. 


-J a. 


a< 


< 




W(\J 




u 


X 


<r- 


w 




CM 


XX 


Qoo 


K •» 


XX 


X 




K* <-IX 


< 




•» • 


to 




z 




jK • 


^ o 


o 




1 


-J-W- 




X 




xo 


O 


X 


< 




•►O 


►- 


uo 


h- 


X 


LUQ'fVJQ'^ 




X’S. 


m 


u 




K 


X «k 


XCsJ 


Mi «k 


z 


i/) . 










< 


< 


'tx-k IL 




ZO 


a: 




1 ^ 






O •'X 


•►om 


It 


a 


II 


»«. •ooo— m 


——CO 


ou 


N-i 


00O>0— •» H-— '< 


II moo^nj 




•»-«i 




in^»-4 


•>in ► 


OCO-^ 


M4UJ 


u. 


•k • 




<K* 


— njx< H CNjin 




P-4QE: 


—1 


inx 


>t<3C 


H-00 




^in^w:^ II ac •ujx 


•►^XX UJX II o 


•^UJ 




II «k»M 


ujr^fM< 


QC 


UJ 


M 1 


M 


— OrcvJO 


•►>o---x 


sO^O 




— 00 


XO'O'^XXnO— 




OUJ 


X 






•^K>r^'-»-2<x II 


wZ 




^HZmK-O 


ax 


H- 




Lli^ 


'-II II — 


|^UJ< •►^UJ<«— 


a< 


Mr ^ O O Ui Ml UJ < < ^ UJ < i-4 < 


K 




(7'X(nOUJK-^Xl-H-l-XXfVJ»-X>--*-i<X • 


• to H ►— X O ^ H- X ►»> O X ^ 


00 


a: 


w 


wOUJXOZO^Q^OO 


N-4 0C 2 0^ x-^o Zi-- QC X 


«-ia:2<QC,j 


•mQC 


O 


OOO H II X-'<O<0CO'^O0COOO<< II 


II XOQ^omocr:ooujo< 


XO 


LL 


03co<rcoKxaruu-3LL < 


O3u.oooxcr<coau3u- 0 3ru.ua:u.o 


►-LL 




































ro 


o 








in 


tn^ O 














ro 


fSI^ 






in 


o 


^ >r 














OJ 


rsjcvj 

















oou 



uouuuo 



113 



1 



1 



i 

I 



K 

O' 















-0 






iL 
























> 


















< 


rsi 




QC 


















X 


•> 




X 




41* 




rvj 










< 


-5 




o 




rsj 








A 






X 










* 




» 




rg 


— * 




X 


UJ 




o 




# 








» 






N. 


00 








4>» 




n 










♦- 


< 




o 








w 




«>» w 


C\J 




X 


X 




Hi 




•V 




o 


(M 




« 




h- 


O- 




M 




o 




2 




-'X X 


* 




9> 






|M4 




X 




w 




Ka a 


N>4 










«J 




w 




1 




-j-* 


N>4 




-0 


-5 




<I 




1 




«>ik 




OiZ 00 


o 




w 






X 








a: 


a: 


o 


W 




X 


CO 




a: 




QC 




3: 




wOO o 


4* 




a 


or 




o 




z 




* 


V. 




*■>4**. 




i-M 






2 




* 








# 


rg— • 




QC 


#>» 


















rvi 




O 


-0 




UJ 




.J 




w 


w 


%0m. ^ 


•K- w 




00 






X 




w 




o 


o 


r- 






X 


< 




K 




o 




1C 


:s 


w •» ^ 


O-J 




c^ 


X 








z 




« 


* ^ 


>-r^ h- 


wLU 




•> 


< 




00 




* 




• 




oc •» •* 


--O 




X 


X 


— 


X 




t 




nj 


nj ^ 


«— < 


>-♦ 




in 


•- 


O' 


a 




eg 




w 


w W 


00 ^ •“ 


QC-* 




»<*> 




w 


< 




w 




1 


1 X 


r- 


O-H 






-0 


UJ 


ct 




1 






# 




00 W 




-0 


w 


u 


o 




A 




rg 


X o 


fSj^H •— < 


wK 






X 


K 






eg 




♦ 


^ o 


-fr • •“ 


♦ ^ 




OfSJ 


Ql 




a 




* 




♦ 


1- * 


* h- h- 


^UJ 




3: 




I- 


z 




* 






-j Qt: 








X-o 


QC 


N-4 


< • 




QC 


in 


zr 


UJ 2 




-^w-5 




in^ 


O 


•k 


o 


in 


1C 


(\i 




o « 


►— 


•-HV.W 




•.yj 


00 




00>— » 


nj 


V 




« 


1 ►- 


lurvi •►rvir*- 




CD 


xoo 


•» 


X 


UJZ 




» 


• 


-^<-11*. 


^ a 


CO ^ •* 


l-LLX 




Q0< 




w 


1-0 


• 


<Wk4i» 


o 


-^Osi 


* 




UJCL< 


< 


•^X 


-5 


UJ 


XX 


o 


—eg 


♦ 


cD«^ X 


X •-r- 


XXX 


w 


-50. 






0.0^ 


* 


x» 






H- w 


**^1^ 


|-<V 


rsj 




QCO 


H- 


x< 








►-ci:<-»— •>. 


<( •> •> 




2 




3: • 


|M4 


ox 




Ko: 


1 






X ^ ^ 


00* "> 


< 




►O 


H- 


o 


1 






UJw 


rgonjo 


w #»r^ •> » 


O 


K 


X •» 


•OfM 


|M4 


o 




UJv 




o 


2w*< 


^in ••inr*^ 


oi^co 


< 


^X^ U- 




zx 




o 


Q04»>C 


I 1 # <W# # 


•«-» 9*^ 9 •* 


w 90. 


II 


••Oao-^in 




0-* 


X*s0 1 * 


•o • 






II o ^o«-< 


II o-^ 


— * 


m 1 -^ ^ 


•*in » 


OX 


— X 


•‘O 






<H-xrsi 






rg 


r^x ^^co'O 


in< 


>->- 




l<— |4i« 


II in II 


wX II 


Lu->< •‘<1^01 n o 


*UJ 


»»»g i>» 


II •« 


LUh*fM 


QC 


II in II 




-> 1 






-^3>O'^-9-0vOw 


X ^ 


om 


.J 1 XK-^ 


II 


-Jh-Xi/5 J II 2<N< 


CO-0 II 






'-'H-zinh- 


Q.X 


II 


-IK 


rn-^O'tU^ 


W IIOJ I 


NMwxniX - 


II 


>-#UJi-#UJ<aDC' LLK 


N-IW< 


h- 




^•^inOUJ>--»QXh-l-i>0 •‘OOCM 




,i-i#t/)KKXO< 






00 


rvi.j(>JOiu 




wOOiX^OZU II ^ II •‘I 


-iQCXwOo<2*-<QCI 


►^QCZ<QC 


*-NQC 




4^0 


OOO II II X— <0<-»^0^>-'U<CD II 


II xoa:oini 


□ OCOOUJO 


xo 


OOO H II 


03r0<CDI-X XOU-O o < 


lD00XCf<®Q.O3CU- ( 


Q3LLOCru. 


h-LL 


03040 






























1 






m 


in 


ino' o 












in 






h- 


00 ^ in 

































oooouu 



















t 


-J 




























! 


o 


























m 




•k 


























•k 




CO 




1 






















— 








' 
















— 




•> 


— 




o 
























«— 






•k 


























UJ 




o 


Uj 


















fo 




w 


00 




•> 


> 




— «k 


















< 


< 




o 


QC 




CM 














h" 




X 


X 






z> 




* 










fO 




o 




< 


Q. 




o 


o 




* 










-w- 




UJ 




X 












— 














— o 




X 


— 




o 


C3 




-0 














— * 




h" 


— 










-0 










f— < 


w 


— »—k 




•» 


— 




o 


o 






— 










l-M 


0sJ(\i 




X 


— F 




•k 


UJ 




o 


OsJ 




C\J 






X 






h" 


o 




o 


KJ 




3 


* 








-IX 


O. 


*h- 






o 




•k 


*p>4 






* 




# 




LUQ. 


w 


1-4 J 






•k 




LU 






1 


— 








Qw 


00 


*-*UJ 




•J 


— k 




•J 


< 






— 








--Z 


O 


xo 




W 


— 




h- 


X 




QC 


QC 










o 


W* 




O ’ 


»— 4 




— 


oc 




3 


3 




3: 




—00 




4- — 




Q. 






h- 


o 




* 










— 




— OsJ 






< 




— 


z 




— k 


— *k 




— »■ 




• 4h 




OsJ* 




0^ 


X 
















o 




O 


— » 


♦ * 




o 


< 






UJ 




“5 


"0 








w— k 


h" 


* — 




00 


X 




UJ 


X 




w 


w 




Q 








N-4r-4 




X 






CO 


h" 




O 


CD 




X 






h" 


Tw 




cr 


— 




< 






3 


3 








Oh" 


•» 


*-'h- 




•k — 


— 






OO 




* 


* 


— 


• 


««» 




IF^ 


-^-1 




XV, 


M 


— k 


•k 


X 




• 


• 


— k 


CM 


f-4 


# •— • 


•» 


KUJ 




tr\ •* 


W 


O' 


o 


a 




OM 


OsJ 






w 


^ »k 


h" 


o:Q 




•N«i^ 


u 


•— F 




< 




w 






1 


X 


fOh" 




OWF 




— m 


a 


UJ 


o 


oc 






1 


X 


«-*» 


■w* 






00 w 




-j •* 


— 




•k 


o 






A 


* 




o 




•» 


wN^ 




^-j 


OL 


1— 


— 






fsJ 


X 


o 




o 


— •' 


ch 


* • 




o^ 


u 


— 


fO 


o 




* 




o 


K 


* 




•» 


— •-4 




XUJ 


00 


h- 


•> 


z • 


if\ 


* 




* 


«J 


QC 


h- •* 








Xoo 




— 


—1 


<u 


OsJ 


OC 




QC 


LL' 


X 


UJ^— -J 




Lh< 


— 


»■ 


w 


1—4 


— i 


3 


UJ 


3 


o 


* 


X •— ^f\j— .1^— .W 


«— 


•-X 


1— 1 


«— 


UJ 


ooz 


• 


— F 


o 


* 


1 


1—4 


*f\jh- 


^w0sJ<I 


OD 


XQ- 


1—4 


00 


00 


LUO 


o 


* 


1 


— 




a 


X 


•*00 


•*!-* X 


•» 


oox 


w 


«— F 


< 


KX 


# 


F— k 




0. 






OfVih- ►h"LU* < 


< 


•<D 


o 


UJ 


X 


:dcc 


— 


— 


—4 


* 




X 


<t •* 




►XU.X 


— 


-J— « 






a 


Q-< 


— ( 


X — 


COw 


X 






X^h" •‘h"h"Q-X. 


OsJ 


X 


^00 


h- 


•k 


XX 


1 


wOsJ 






<-J— 


X, 




•*CNi 




z 


f-HX 


1-4 • 


— 


o 


a 


-5 


h"* <-J— 


V 






wCNIiM ^pH00<O 


< 


•*o 


— o 




X 


ux 


-> 




wUUfvj 


«4fc 


OsJCOsJO 


w ^ 


•*m 




h" 


X<-4^ OsJ 


— 




h" 




LUQCOsJO* O 


< 


mh^ *h-o* o 


< 


>4“ •‘00— u. 




y-4 


ZQC 


GO 4- 


Q3Zw* < 


* 


II • 


•*o 


•*wira. 


II 


— 


•*irsm 


— — 


kOO 


03 


1 ^ 


<-F — 


* 




• 


II 


— 




OCO 


I w 


— O 


m4 • 


•*— * 


►— 'X 


• 








•*— Co: 


ro 


a^x-- 


II 


s0<X 


KtL 


II lf\«— 4i— 4 f«— 


<[h-^ 


OsJ 


HOC'-' 


— 'Uj-J<roo(<>-J II o 


•*UJ 


I •>F-i40* 


— * ^1-4 LUh" OsJ <t 


QC 


3 1 


II wX II OC^ 


k— >1 


—OK 




D •»-! 


•*< 


•*^— t/) 


-iz?v0--e* 


-4 0^: 


O •-•0^ 


OUJ 


3 II Xh-w 


^o-j 1 : 



II II 






w. n Lu i-Mwx •“V) «-»uj*-^uj<ir\h-UJ<‘-^^< 

h— *QXl-h-*/>fVJ II OsJ»-< j^«-i</)h-h-z •‘>O^XKOZ*J 

Lu:^wuzuii •‘o:xwxx<z*-‘QCx ^ccz<ioc^ 

II M xoa:oa'OQCOOLuo< 
t-x XUU-X X coXO<COQ.U3rU- Q3:U-OQCU..U 



CLX 

I/) 

^Qc: 

xo 

h-lL 



rn 

r\j 



f- 4 i/> 

OsJ O' 



vr\^- 

0^0 



o 

>0 



*X(/)UJ H Z 
f^-XMJJ-Jw H Q ^ 
fM"50sJOtUh-— wXK 
^ ^OUUX OZ 
OOO II II X'-'-f <o 
03CQ<CQh-X XO 

O' 

fM 



UUUUUU 



115 



















00 














h* 










► 












•> 


•> 


* 




•» 




'O 












r- 








—* 




► 








T 




•» 




m 




X 




0 


UJ 




, 


X 




f-^ 








w 




► 


> 




A 








•> 


1- 




< 




0 


ex 




rg 






r- 


r- 


-j 




X 


X 




x> 






X 








UJ 




< 


w 


0 


0 




« 








«h 4 






X 


UJ 


► 






A 






•» 




— * 




X 


00 


0 


Of 




r> 














f^ 


< 








w 










rgrg 






X 


0 


0 




W 






rg 


'-^rg 


4 h^ 




X 


a 


► 


UJ 




0 






W ^ 


M » 


♦ 1- 




r- 




0 


M 




a: 


K 




•-irg 


xrg 


•—-J 








► 


w 






UJ 




X •» 


a 


*^ilJ 






X 


UJ 






1 


X 




am 


'--'m 


cxo 




X 


w 




< 






H 




»«-» » 


00 •* 






w 


0 


K 


X 




(X 


X 




z-^ 


0-* 


4. 4 P-» 




0 


a 




ex 






r** 




NH ^ 


o 


— m 




a 




K 


0 




* 






com 


m 


rg-tt- 




w 




w 


z 




A 


X 




•» 


* •► 


« ♦ 




ex 


X 








"3 




X 


♦ iTk 


irv 


* — 




0 


w 




UJ 








o 


-—rg 


-^rg 


^•-4 




00 


< 


UJ 


X 




W 




< 




fWf-^ 


q:^ 




X 


X 


CO 


K 




0 


CO 


X 


► • 


^ • 


wK 




(> 


< 


< 






X 


X 












X 




00 




« 


f-^ 


A 




KUJ 




XV* 


•» 




X 




• 




X 




QCO 




in •* 


«-* 


cr 0 


a 




rg 


X 




•‘-J 




Ow 






X 


w ► 


< 




w 


00 


X 




00 ^ 




— >t 


w 


UJ 0 


ex 




1 


«m 4 




►X 


•‘X 






X ► 


• 0 


-j ► 


0 






•» 


««» 


^cor-^coh* 


«• • 




•^x 


a 


K 






rg 


< 


X 


•» 






-0 


Ow 


1—4 


— ^ 


Q 




« 


X 




m 4 ^m ^ ^ 


— . '-i' 


-5 


3 :uj 


ex 


0 - ► 


z 


in 


# 




h- 




^ <► 




w 


XOO 


0 




< 


rg 


ex 




UJ 




*** 




< 


— in< 


00 


► w 


0 




3 C 


X 


X 


»>*— » 


•» «M>i» •» 


%--m 


X 


CO •‘X 




►► UJ 


00 N- 


• 




o 


H- 


fsir^^cgr-^ 




< 


xa 


A 


00 00 


UJZ 


0 


■H' 


•«4 


•k 


■w ^ 


^ 


UJ^f- 


X 


< cox 


X 


w < 


1-0 




A 




CD 




XU. 


V 


W *0 


w 


UJ X 


OX 






X 








c-»a 


A. 


rg x^ 


0 -- 


«j a 


aa 




X— 


X 


<f^ 




---X 


-5 


Z X 


3:00 


h- ► 


x< 


1 


^►^rg 


*>-< 


^ • 


^ <► 




00 < 


-0 


< ^x 


•» • 


NM Q 


ox 


-> 






xo 


m^rgm.^rg 


0 -- 


w 


K •‘O 


xo 


►- :s 


0 


^>4 




X'^ 


rg 


•» ^ 


» ^ » 


0* 


0 


< x^ 


rg 


NH ► 


X 


f-^ w 


uuex 




m>togin^rcM 


win 


a 


II 


Fl-^U. 




Zl- 


CX)^ 


ox 


O •‘X 


►irUTi 


• «M| 




II • 


w 


^ in 


oDinm 




ou. 


-O'O 1 — 


-tO{j^>t •* 




^0— 


>QC 


^ ^^x 


►rg ► 


OQO— » 




t 




rox< H rgin 


— ►X— ►T 


-> M — 


•0 


•» ^ X 


p^^g 3 


r*< 3 : 


l-U. 


II 




•00 


->LU •‘^011 •H-<ujr^rM< 


OL 


^ 1 


II --X 


^o-^x 


vO— 3 _l •• 


►'-'QC 


: — )>t II 


-^oowax>o ^3 


► #-^eX 


OUJ 


*- II XK<^ 








^Z-^l-X 


'-^^-zini-Q 


ax 




-Jl- 


UJ< 




•X »#«x ^o If a^H -*5 


«ujM-uj < in UJ < ►-I 


1W< 


H- 


m'^r*uu J 


^ X X ^^J K X h- </^ fv» 00 1/) fvj o 




•^OOI-H-X ►mh-XH-OZ-l 


oo 


mi-imouj 


mOCOO 


•--QCZH 


M ^ 


• CXXUJwqcoKZ'-QCX 


>-exz:<cx-i 


►hCX 


w 


wQ 


a:o-' 


OQ^OO 




' 0<00 II 


II XOa:oa'Oa:oouJO< 


XO 


000 M II 


3 U. 


osu-ooe: 


(X QCOOX 


Of<coaoxu, 1 


O3u.oau.o 


l-UL 


O3CO<c0 






^rg -«rg 


















o 










r- in 


in^H 


1 0 








sT 










eg f-g 


rgm 


h* 








CM 


rgrs» 












h* 






f 



oooooo 



116 



o 



(VJ 


A ^ 








00 












h* 








•> 










« 


•» » 




















f^ •* r^ 








•> 










:ac 


•«fVJ •. 








o 


lU 






» 


w 


1 

♦ 

1 

t 




in 




► 


> 








H 


►rg •‘Cg 


< 






o 


QC 




eg 




U 


ir> •* if\ 


X 


NM 




•> 


Z5 




# 




lU 


►fM 


< 


<1^ 




o 


O 




» 




O 




X 


uu 




p> 






i—i 






•►-4 •►st — * 


X 


CO 




o 


car 




X 




+ 




r- 


< 




H- 






X 






►'4’ ^ fvjfsi 




X 




o 


O 








fVJ 


^c\j » i-^rg ^ ^ ^ 


X 


a 




*• 


uu 




Q 


rg 




^ •‘U> — X 


h- 


•> 




o 


rvj 




3 


*w^ 


♦ 


Xfvjrg— *Q.rg j 


•> 






•> 


»»4 




«>»» 


* 




a. » •h-'^ ►cor^*— fLu 


iiiiik 






uu 


«J 






— 




►</>>t(M •‘1/50 


k-4 


•«pi> 






< 






-** 


o: 


X •h-^* 




UJ 




h- 


X 




QC 


X 




•—••—< <r »o*^o 


uu 


a 




hh 


QC 




X 


:r 




CO -^'4‘ 


a 






h- 


o 










>t< 


pmt 






1— * 


z 






— • 






QC 












X 


X 




iT^co •— •-* — 


U 






•J 


IXI 




X 


X 




-^rg mh^rgxu^*— ^ 


</) 


< 




UJ 


X 




w 


w 


Q 


^ 1—4 {/) ^(y) 


X 


X 




X 


♦— 




Q 


o 




^ tlTi •^^wH- 


o- 


< 




< 






-X 


X 




iTiO*— “ •* ••O-^ •‘^—1 




X 




•J 


CO 




# 






#. w -i^ fc> * fO H- CU 


XV^ 




««k 


» 


X 




• 


m 


rg^ 


^<h- •► **<r^ •*a:c 


IC\ H 




0^ 


o 


o. 




rg 


rg 






»M«i»» 


pmt 






< 




•w 




1 X 




-^u> 


w 


UJ 


o 


QC 




1 


1 


— * 


•»X ••rg^x ►rg^v^— ‘ 




UJ 


u 




o 






Mi 




•. •►COh* •-*.«— 


■■Kl 1 


o. 


H- 








rg 


X 


wO 


•^fO •►rg-^^-i-.*^ 




HH 


HH 


in 


o 






w 


K* 


rT)4*i— ». 4» 


2UJ 


o: 


K 


•> 


z 


in 


*»► 


1- 


UiOC 


p, 


XCO 


u 


»-* 




< • 


rg 


QC 




UU3: 


‘K< 


— \rs< 


CO 




w 


u 




X 


UJ 


o* 


►arg— ►cLw^t^jx 


CD •‘X 






UJ 


COh-i 


• 


w 


o 


1 M 


rgi^,-4'-* ^ uj< 


•* xo. 




QC 


CO 


UJ2 


o 




1 


-‘O. 


•» •* •-Z4' •‘OOUj-frOX 


< cox 


»»4 


w 


< 


♦-0 












w hO 


w 


UJ 


X 


IDX 








ao-^x 


► •- •</)»— 1 •> ^Oh-Q. 


rg 


Q-. 


«J 


Ql 


qlo: 


»— 4 


X— ‘ 


CD^ 




••tn^ wX— 


z: X H 


xoo 


K 




x<r 


1 


•^rg 


►1— 


<-i>^ 


► »• >J' ^ (/><[ 4" 


< r-*X 


•» • 


•— 1 


o 


QX 




K* <l_l 


wlL‘«-^ 


>^,^rg ^i^rg Qw—w 


h- »0 




K 


3: 


u 




•J» 


•-CU 


rgOO 


». p, •» 


< X— » 


rg 


»-4 


•> 


X 




ujorrgQ 




ir»>trgh- •’4'rgr^'-'in,ja 


1) —’4“ 






Zh- 




axz^ 




• ►O ► II *10^ 


— \t\ «i— *< 


COUfMTS 




OD 


ox 


OO 1 <- 


Cw 


• 


o ^ Hi* r*- •«*-*-* H-# h- o o 0^ 


UT\ CO^*^ 




oco 




»-i4»«4 


!l • 


1-^ 




^ •‘X •►< •‘X H-0 11 * o 






Q0<X 


HCO 




<t h- 


II a:'^uj<fncv^— 


--)iij ►^uj II ••HHUjr^rg< 


cc 


X 1 


II '--X II oc 


-^O 3*J ^ ►w •►w-^ro II 




N-vo^r5 •—* 


^QC 


aui 


X M XI— ^ 


-^car 


^</) II z<^zir\a:^(/)iryuj*-iw— 


^Z^h-X 


—M-ZiCh-O 


O.X 


A 


-IKXCO 


M ►-•X »C/> ►O ►CL^h-"^ 


UJ H-I UK O' uu < 


iw< 




H- 


^xmuj*j 


^ II 


Xh- corgcoi-4 ii cvjo^*-'<*j^h^ 


•—ooK^-X ••( 




CO 


^■^'i'OUJh- — 


ULJXO^II •» ►o:XUJ— 'COOO<Xi«iOCK 


►->a:z<a:-j 


H-IQC 


«Wi> 


wQLUi^ 


X^<0*-I^* fOO<QUJ II 


II XOQ£O0'O0COOUJO< 


xo 


OOO II II x^ 


►-XXU00 CO cox ar<CDCLU3:U- ( 


D3CU.OQtlLO 


I- CL 


QXQ<XHX 




^rgrn ^ogr^ •-« 


















h- 




fO IT> 


in O' 


o 


\ 










m 


m fo 




CO 


















h“ 













oouooo 



117 



(DELT(K)** 2 H 



o 



CD 



X 

# 

o 

o 

♦ 

QC 

* 






OiC\^ 



inm— ► 

•w- 

^iTv •» ►rvj » *-<irv •»rvi 

^ • *1 

w •o^O'-'cnor^X •►oo^nof^cNJCNJ 

H-W>r\Jw •*1— ^Q^fVJW W 

^1-- •*<«-^XrsJ<^*»- 

UJX<NJ^ •►a •*-! •‘oofvj^ •‘a ••^ 

QQ- •^<in^in<oo ►<trN^ir<vOi-iijj 
wwif\X •‘X •‘X •►OlT'X *“0^ -“X ►KO 
V>X •^^*-1^00^ ►C0*-^0'^00^*—’* 

— •*</^ ► ^ •‘O •“ •‘ 4 '—* 

—00 •>^rT) iri 4 *fO •--fro 
• ir» ••«)<“ •• •‘fsj* 

Oit* 4n-^f\js0 in— (M*# 

(\j in^ •►— rsj— » •.in-^ •.^f — 

— rsjf^ (Njf^ >0 rsj incNj rvj rsj ^ 

wvO^ ^ •* ^O ^ •“ •» »■ »o ^ — 

h- •► •p*>»(Ni>o #r^(\j^ •^-^^J^o 
Q^mo ►o •‘O »o •‘^'-'^>0 

»X ^X ►< •“X •*< ^Xct:OK 

^woOwUU 

>c •‘Xinz^XvOzmxmoo^XvOoowNsO 

# ^l/) ►OO »>N-i ».l/) •‘O* 

# •• 1— ioom ^ 00 ^ 

<-*P 0 4 ^ **4 ► •*4 ► ^4 •“ 

< •■ ^»(\j >t*rsj >t*fvj #ww 



•rsj^^— m 



UJ 

CDCNjf^ ►h- ••h- 
* ► 

Xm>o •►>0 •‘P^ •►N' 



rsj-— win-jx 
..fs. •■r-h- *uj< 

(\i ••in ••ro •-uj*ox 

•»s 0 ••'O ••f^ •‘f^XU.^'v. 



O ► MNJ ►m •*(%] 

*v X •> •‘in ►*-< ►m ►oo^w-^ 

— win^x^r^ ►^X'O •—ixtn •*^x>oo^>--' 

O ^ ►Q- •‘in •‘CL *0 •‘CL ••IT ••CL ••O* -JUL 

< wiri>tw^(M>t^<Mw^w^rsj>tw.(MwinujCL 
♦ H • ^Z •■ • ►Z •</) •• • ^i/> ► II •O*-^ 

• ‘O* oc 

CNJ ^</) ^ 0 .) ^O M O 

wuji^<fn nj<rn rvixm <\j<fn rg^-- ^oo 

3 •‘-J *4^ ••-! ^* ^40 •‘^ •‘-j ••# ••— It 
H ZvO<^ rsJ<#M fNj ^ (M<^ 

••X •— •‘II •— ►x •‘— •‘CLwh-X 









>0 


— 






•k 


N -4 


— 




0 








» 


< 


•k 




0 


X 


X 




•k 


< 


•— >■ 




0 


X 


UJ 




•k 


X 


00 




0 


p^- 


< 




•k 


•» 


X 




0 


X 


a 




» 


r- 


•k 




0 


•k 


— » 




»■ 


— 


X 




UJ 


Ni 4 






-J 


W 


U. 




K 


u. 


or 




M 


CL 


•k 




►- 




— 




N - 1 


QC 


X 




» 


0 


W 






00 


< 




UJ 


X 


X 




CO 


O' 


< 




< 




X 




-J 


X>K 


•k 


— 


•k 


in •* 




cr 


0 


•k— 


X 








w 


uu 


0 




u. 


-J 


» 




CL 


K 




Ow 




»-4 


nO 


ZUJ 


QC 


K 




xoo 


0 






in< 


00 


•k 




••X 


•k 


— 


UJ 


XCL 


— 


CO 


oO 


COX 


X 


W 


< 


•‘O 


w 


UJ 


X 




0-^ 


-J 


a 


X •* 


3:00 


K 


•k 


- 4 X 


» • 


l-H 


0 


•*0 


xo 


K 


3 C 


XrH 


(NJ 







fVJ 

z 

< 

K 

< 

It — •.^— u, ^ 

— m •^ooinm 
vO in^^^ ^>0 ^ ooow 
» ^ ^ 3 K 

xuj ••fHu. II •••«ujr-(M< 
igCOO'-'OrXxO'-'X •‘•-QC 
wZwi-x wKZini-o 
•- UJ •-• Lu < m UJ < ^ < 



Xi-i-oorvjp-inoofvir-in rvir^inoorvjfwtn*-i<«j^^-»i«ooKKX ^^i-i-Xh-ox-ja 

OZOII •»•••“ •‘ •• •»-• •• •‘ •• •• •• ••QCXUJ'^H-l^<Z*-«a:x *-iqcz<qc«joo 

4<0<*“'^'0^4<Mr-^^-— ^'0-^4<MpN-^U<QlL II « XOQCOC^OQCOO^JO<^-Z 
XOLLI-- ^ OOX O<C0CXO3SU- O^U-OOCILULOUJ 

^ ^fvjm^inoh- -^f\jrONj-in>op^ ^ - 



m 



^ in 
in 



inr^ 

>0^ 



o 

c* 



118 



*md 










♦ — 






















»• 








X 








z 


#x 


tfXk 




«» 


« 


1—4 






— * 


fmmt 


— 


•J 


z 


•> 




•> 




z 




v: 


X-!> 










► 










LU^ 




^x 


M4 


Xk- 


n 






KLU 




*♦ 




X 


M 


* 




♦ h- 












:^Z 


-> 


Z — 


z 


^z 




+ Z 


•k 


X* 


•k 




•J 




z— 




h- 




•f^ 




^ W 






h-UJ 


X 




•X 


U4X 


->u- 




^ •►U- 


XI- 


•CL— ‘ 


--X 


•— ^Q- 


l-'f 


QZf^ 


-w- <— 


->^-z 




w<*-. 


’-5 




iC-lX 


< H- 


« » 




4-h-iQl — 


-JXUJ 


» 


wUJX 


— ^ m 


< •X 


u.-- 


*-^x 




Zl-h- 


aic 


XL-KOwUJO • 


LOLU 


Z:^ 


CL UJ 


• h-XOO 


XZ 


<♦ 


ZXOUJK>OfO 


Zh-O 


** 


ZOI- II X •►IL 


a 


•O— 





II nz 
h-OZ M <-50^ 

OZtU-^ JwO 

zzzr^<xi-o 
30 *-iwXwUJZ 
U.OQX</> OCLU 



H-n/)zr^* *^KZ 
HZO^^!^iii< 0 ^ 

M 

ZXXOJ^ *-iQCh-p 
0 *-^ 0 XX 4 »QC 0 UJZ 
ULOUH-Ql IXlLOCai 



o 

o 

n 



119 



B-7 : Program FQUANAL 



This program is an adaptation of work done by Roy M. Johnson, 
Assistant Professor, Naval Postgraduate School. Basically the program 
takes a waveform obtained directly from the cavity and performs a 
Fourier analysis on it. The output of interest is a set of Fourier 
coefficients and their associated phase angles. The program is 
written in such a manner that a table of definitions is not appropriate. 
Each variable is defined within the program. 
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Pressure Waveforms Inset is experimental result, — -sine wave, 
FINAMPI prediction, FINAMPIZ prediction 
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Pressure waveforms Inset is experimental result, sine wave 

FINAMPI prediction, FINAMPIZ prediction 
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Pressure Wavef oms . Inset is experimental result, 

sine wave, FINAMPI prediction, 

FINAMPIZ prediction . 




Q- curves for the fundamental, M = 0.004, 

QCURC prediction, QCURVES prediction 

Experimental results are indicated by ## 

Figure C.5 
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Q-curvesfor the second harmonic, M = 0.004 

QCURC prediction, ----QCURVES prediction 

Experimental results are indicated by 

Figure C . 6 
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Q-curves for the third harmonic, M = 0.004, • 

QCURC prediction, QCURVES prediction 

Experimental results are indicated by •• 

Figure C.7 
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Q- curves for the fourth harmonic, M = 0.004, 

QCURC prediction, QCURVES prediction 

Experimental results are indicated by •# 

Figure C.8 
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Q-curves for the fifth harmonic, M = 0.004, 

QCURC and QCURVES prediction 

Experimental results are indicated by •« 



Figure C.9 
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QCURC and QCURVES prediction 

Experimental results are indicated by#« 

Figure C.IO 
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Q- curves for the fundamental, M 0.005, • • 

^QCURC prediction, QCURVES prediction 

Experimental results are indicated by#« 

Figure C.ll. 
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Q-curves for t^ie second harmonic; M = 0.005% 

QCURC prediction, — .--QCURVES prediction 

Experimental results are indicated^ by •• 

Figure C.12. 
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Q-curves for the third harmonic, M = 0.005,^ 
QCURC predict ion ,j;^'-QCURVES prediction 
Experimental results are indicated by •• 
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Figure C.13 
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Q- curves for the fourth harmonic, M = 0.005^ 

QCURC prediction, QCURVES prediction 

Experimental results are indicated by •• 
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Figure C . 14 



136 




Q- curve for the fifth harmonic, M f 0.005, 

QCURC and QCURVES prediction 

Experimental results are indicated by#«. 

Figure C.15 
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Q-curve for the sixth harmonic, M = 0.005, , 

QCURC and QCURVES prediction 

Experimental results are indicated by## 
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Figure C, 16 
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Q- curves for the fundamental, M = 0.009 

QCURC predict ion, ---QCURVES prediction 
Experimental results are indicated by#. 
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Figure C.17 
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Q-curves for the second harmonic, M = 0.009, 

QCURC prejiiction,----QCURVES predi'cticm 

Experimental results are indicated by#. 

Figure C.18 
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Q-curves for the third harmonic, M = 0.009, 

QCURC prediction, QCURVES prediction. 

Experimental results are indicated by#. 

Figure C.19 
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Q-curves for the fourth harmonic, M t 0.009, 

QCURC prediction, QCURVES prediction 

Experimental results are indicated by #« 

Figure C.20 
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Q-curve for the fifth harmonic, M = 0.009, 

^QCURC and QCURVES prediction 

Experimental results are indicated by 

Figure C.21 
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Q- curve for the sixth harmonic, M = 0.009, 

QCURC and QCURVES prediction 

Experimental results are indicated 

Figure C.22 
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Phase dependence of the various harmonics on^4P 
as predicted by PHAMPC. n = harmonic number. • 

= data from FOUANAL for n^-th harmonic 

Fi,, '' ->3 



145 



INITLA.L DISTRIBUTION LIST 



No. Copies 

1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 22314 

2 . Library ' 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Commander 1 

Naval Ships Systems Command Headquarters 

Washington, D, C. 20360 

4. Commander 1 

Naval Ordnance Systems Command Headquarters 
Washington, D. C. 20360 

5. Professor Alan B. Coppens 8 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 

6. Professor James V. Sanders 1 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 

7. CDR Wayne "L" Beech 1 

1108 Leahy Road 

Monterey, California 93940 

8. LT Paul G. Ruff III 1 

USS Dale (DLG-19) 

FPO San Francisco, California 



146 



UNCLASSIFIED 

Security Classification 



DOCUMENT CONTROL DATA - R&D 

(Security ctmeeHication of title, body of abstract and indexing annotation must be entered when the overatt report is classified) 


1. ORIGINATIN G ACTIVITY (Corporate author) 

Naval Postgraduate School 


2a. REPORT SECURITY CLASSIFICATION 

Unclassified 


Monterey, California 93940 


2 6. GROUP 


3. REPORT TITLE 




A Theoretical Investigation of Finite Amplitude Standing Waves in Rigid Walled 
Cavities 


4. DESCRIPTIVE NOTES (Type of report and inctusive cfa^e«^ 




Master *s Thesis, December 1967 




5. AUTHORC^^ (Last name, first name, initial) 




RUFF, Paul G. , III 





6. REPORT DATE 

December 1967 



8«. CONTRACT OR GRANT NO. 



b. PROJECT NO. 



d. 



7m. TOTAL NO. OF PAGES 

i44 



7b. NO. OF REFS 

26 



9a. ORIGINATOR'S REPORT NUMBER('S; 



96. OTHER REPORT NO^S^ (Any othmrnumb^rm tfiat may 6a m9»ignmd 

mim rmporO 



10. AVAILABILITY/LIMITATION NOTICES 



11. SUPPLEMENTARY NOTES 



12. SPONSORING MILITARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 9394^ 



13. ABSTRACT 

The Coppens- Sanders perturbation solution for the one-dimens ional non- 
linear acoustic wave equation with dissipative term describing the viscous 
and thermal energy losses encountered in a rigid walled, closed tube with 
large length- to-diameter ratio was extended to include sixth order terms. 

The solution was then investigated to determine the region of validity. 
Computer programs were written to evaluate and graph the resulting waveforms. 
Available experimental results were compared with the theoretical predictions 
and good correlation was found to exist in the region of low Mach numbers. 
This agreement was found to gradually deteriorate as the Mach number was 
increased. A Fourier synthesis approach is also presented and the leading 
terms of the first ten harmonics are derived. 



DD 



FORM 

1 JAN 04 



1473 



147 



UNCLASSIFIED 



Security Classification 



UNCLASSIFIED 

'' Security Classification 



KEY wo RDS 



ROLE W T 



LINK C 



Finite- Amplitude 
Standing Waves 
Perturbation Approach 
Fourier Synthesis Approach 
Rigid-Walled Cavity 



DD ,'°“..1473 'BACK) 



S/N 0101-807-6821 



148 



UNCLASSIFIED 



Security Classification 



A- 3 1 409 



